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Within the framework of a modified Balitsky-Kovchegov equation, we calculate 
and provide estimates of non-linear saturation effects expected in the LHC range of 
energies. 

The main objective of this letter is to provide reliable estimates for high density QCD (saturation 
effects at the LHC range of energies. Our estimates are based on the ^alitsky-Kovchegov 
equation which is the mean field approximation to high density QCD 

ail. 

In spite of the 



restricted theoretical accuracy of this equation 0, Iz, 0, S ll^ llll. 13, 3, 15, 16, [l7| in 

the framework of the colour dipole approach the physics underlying it is rather transparent, 
and it has led to a successful description of the main features of the experimental data, both 
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for deep inelastic scattering and ion-ion collisions 

the Balitsky-Kovchegov equation is the theoretical understanding of its mathematical structure 



Among the advantages of 



the analytical solutions in the restricted domains |2 5L l2 6 L l27l |. as well as extensive 



experience in numerical solutions 
moment the best that we have. 
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This equation is at the 



However, the Balitsky-Kovchegov equation has very serious deficiencies: it does not include the 



next-to-leading corrections to the BFKL equation 



39^ which are large and those should be 
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taken into account to obtain reliable predictions at higher energies j41|. In the estimates given 
here we rely on the modified version of the Balitsky-Kovchegov equation which was suggested in 
Ref. 4^. The modified Balitsky-Kovchegov equation includes the correct re-summation of the 



NLO corrections to the BFKL equation suggested by the Florence group 43]. It also contains 



the important observation of the Durham group 4^ , that the main part of the NLO correction is 
the inclusion of the leading order anomalous dimension of the DGLAP equation this in the 
framework of the BFKL approach. 



The modified Balitsky-Kovchegov equation suggested 
dN{r,Y;b) Cpas f d'^r' 
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En. (HI) includes the full anomalous dimension of the DGLAP equation in leading order, in the 
approximate form proposed in Ref. |4a|, namely, 



as 
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which describes the exact 'j{uj) within an accuracy < 5%. The advantage of Eq. is that it 
conserves energy even when including the non-linear term. 

Unintegrated structure function: 
Solving Eq. (P) we obtain the dipole scattering amplitude. The more transparent physical meaning 
is the so called unintegrated structure function, which is the probability for a hadron to have a 
gluon with fixed transverse momentum. This function enters the calculation of the inclusive 
production for the gluon jet, as well as all other inclusive cross sections. 3, 
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dy d^pt — ^ Pt 

where the unintegrated structure function (p is determined from the solution of Eq. (P) by the 



following equation 



(p{Y,kt 



N, 



4:asikf)7i 

where is a two dimensional Laplacian and 

N{Y;kt,b) = 
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and the gluon structure function is equal to 



as{Q^)xG{x,Q' 



dklas{k1) <P{Y,kt) 



(4) 
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Suppression due to the saturation effects: 
In Fig. [T] we plot the ratio 

D {kt; x) 



(7) 



4>^ {kt,x) 

where (p^^ is the unintegrated structure function obtained using all terms of Eq.(l) , while (f)^ is, 
the unintegrated structure function obtained by excluding the non-linear part of Eq.(l). 
Our initial condition for solving Eq. are the dipole-proton amplitude from the CTEQ param- 
eterization for the gluon structure function. We use the following form of this initial distribution 



(r, x; b) 



exp 



asTT" 2 
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(8) 



where S{b) is a dipole-like profile, R is the size of the proton, we took = 10 GeV 



500 




FIG. 1: The ratio given by Eq. ^ which shows the effect of the non-linear terms on the prediction for 

the dipole scattering amplitude at the LHC energies. 



From Fig. ^one can see that the non-linear term in Eq. is important, and in the LHC range 
of energy it suppresses the value of by 30%. It should be stressed that the inclusive cross 
section is proportional to 0^, and we expect even a larger suppression for the inclusive production 
(approximately twice as large). Note that kt) < 1 for kl < Q^. For k^ ?a Q^, on the other 



hand, we observe an antishadowing effect on which we shall comment 



be 



of the ratio was expected theoretically (see for example Refs. 
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Such behaviour 



ow j^. 

, however there are 
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no reliable estimates in the approach which is based on the solution of the Balitsky-Kovchegov 
equation, and which describes all data at lower energies. 

Increase due to the saturation effects (antishadowing) : 
Fig. 121 shows the value of xG{x^Q'^) for the non-linear equation (see Eq. (Q)), together with the 
value for the solution of the linear equation (Eq.(l) without the non-linear terms). Note that 
the solution to the non-linear equation leads to a larger amplitude for large values of x. This is 
not expected as the non-linear term has a negative sign, and manifests itself as the shadowing 
correction which decreases the value of the amplitude. The increase of the ratio D (see Eq. (0) and 
Fig. El) at low energies is significant, and cannot be explained by errors due to the approximate 
formula for the anomalous dimension of the DGLAP equation, or due to the numerical prodecure. 
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FIG. 2: The value of the gluon structure function as function of x for different values of . 

Non-linear denotes the curve for the gluon structure function calculated using Eq. ^ from 
the solution of Eq. (0j, while linear is the same without the non-linear term. CTEQ denotes 
the curve for the gluon structure function given by CTEQ parameterization. 



The explanation for the ant i- shadowing effect is very simple. At large kt, the non-linear term is 
small and the linear evolution equation is a good approximation to the solution. Decreasing the 
value of kt, the largest changes are induced in the slope of the amplitude (or/and gluon densities) 
with respect to ln(/ct) (see Refs. 0,0,0)- 

The slope for the linear equation is larger than 
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for the non-linear equation. For large values of kt, the unintegrated structure function is equal 
to the solution of the linear equation, while for small values of kt, becomes larger than the 
solution to the linear equation. 

Another way to see the same effect is by comparing derivatives with respect to Y see Fig. |21 For 
the solution to the non-linear equation at large x the derivatives with respect to Y are much 
larger than the same derivative for the solution of the non-linear equation. Therefore, for the 
case of the non-linear equation the linear term, which contains —d/dY, is larger for the non-linear 
equation, leading to an increase in the value of the solution. The unintegrated structure function 
in Eq. (jH) reproduces the same qualitative behaviour as the dipole scattering amplitude. 
All these arguments can be obtained directly from energy- momentum sum rules which hold for 
the modified Balitsky-Kovchegov equation, namely, that the integral over x is independent of 

j i.,^^ (y;Q^) ^ Co,.sm ^ j (Y;Q^) (9) 

In the region kl < Ql{x) (or x < Xs, where Xs is the solution of the equation k^ = Q^Xs)), 
we expect D{x, kt) < 1 (see Fig. 12)) and based on Eq. ((Tj), we anticipate that in some region of 
X > Xs, this ratio should be larger than 1. 

Fig. 121 displays quantitative estimates for these antishadowing effects. 
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FIG. 3: The ratio {4>^^ — for different values of Q"^ as function of y = ln(l/x). 
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The effect is rather large and has been discussed previously in Ref . . In our approach it arises 
naturally 51 1. 

The comparison of the results of the present calculation with those of CTEQ and MRST is not 
straight forward, since we calculated the unintegrated gluon distribution, while CTEQ and MRST 
present results only for the structure functions. We are presently busy with this project and hope 
to be able to present our results in the near future. 

We hope that our estimates of the scale for the nonlinear effects at the LHC energies, will help 
to produce more reliable calculations for the cross sections of 'hard' processes of interest at high 
energies. 
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